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This Letter proposes a rescaled adaptive coupling scheme for the synchronization of spatially 
extended systems. Coupling and synchronization are analyzed from the point view of image filter 
construction. A length rescaling technique is introduced based on dimensional argument of con- 
. . ^ trol process. Control sensors are adaptively selected according to local information. The coupling 

^^ . strength on each sensor is automatically adjusted according to the magnitude of local oscillations. 

^^ ' We demonstrate that the present scheme can efficiently suppress and control spatiotemporal oscilla- 

tions and thus, provide a powerful approach for shock capturing. Both the Navier-Stokes equation 
describing shear layer flows around a jet and Burgers' equation are solved to illustrate the idea. 
Ph; PACS numbers: 05.45. +b, 02.70.-c, 47.40.Nm 

.^ . Synchronization phenomenon is of fundamental importance in telecommunication ||l|, electronic circuits |^, non- 

linear optics H, and chemical and biological systems B. The phenomenon has been studied extensively both by 
t^^ , numerical and experimental means. It is believed that an in-depth study and understanding of synchronization vi^ill 
CN ■ greatly benefit the advancement of science and technology. Different types of synchronization, such as identical SQ , 
generalized H, lag and phase synchronization H, were proposed. Recently, synchronization and control of spatially 
extended systems have received great attention [pU9). For a given system, the degree and rate of synchronization 
depends vitally on the coupling scheme used. A variety of coupling schemes, such as unidirectional coupling, receptor- 
product coupling, weak coupling, strong coupling, global coupling, and local coupling, have been studied. However, 
P5 adaptively coupled schemes and their effect on the rate and degree of synchronization have not been addressed yet. 

|Jh ■ Moreover, very little is reported on synchronization with respect to the understanding of nonlinear hyperbolic con- 
G ■ servation laws, shock capturing and, in general, computational methodology. The latter has had tremendous impact 
to a variety of disciplines in science and engineering. In fact, much of the present understanding on synchronization 
^-^ • was achieved with the aid of numerical computations. 

^ ' The main purpose of this Letter is to introduce rescaled adaptively coupled synchronization schemes and to use them 

OO for shock capturing. A local gradient based coupling scheme is introduced for continuous systems. Nondecreasing 
\l functions and non-increasing functions are designed for oscillation reduction and image edge preservation, respectively. 

We demonstrate that appropriate coupling of two identical dynamical systems can result in a surprisingly novel and 
efficient scheme for shock capturing. The validity and efficiency of this scheme is tested by using an inviscid Burgers' 
^^ ] equation and the Navier-Stokes equation for incompressible flows. 
f~*) ' For simplicity, we consider an identical synchronization, where two coupled systems are exactly of the same type 

"^ ' and can be given by a partial differential equation of the from 

'-^ ■ ut = F{u,Ux,Ua;xr--)+c{u,w), (1) 

^ , where c(u, w) is a dissipative coupling term, which is proportional to the difference between the states of two systems, 
l^ {u — w). Based on an experimental consideration, Junge and Partitz B proposed a sensor coupling scheme, which 

'^ utilizes the difference between two spatially averaged local signals (u — w). Here, w{x,t) — j /^JLj^ '"^(j/jO'^y i^ the 

5r . local average of w over a length I at the position of a sensor. It is noted that a dimensional scaling effect has to be 
accounted when the solution of a partial differential equation is related to an experimental measurement. In fact, 
in control experiments, two coupled systems might have very different spatial extensions (the control part is usually 
much smaller in size, e.g., a human brain vs human body). Mathematically, the dimensionless partial differential 
equation usually has one or more parameters which contain the length scale. We argue that, when two dynamical 
systems characterized by different values of a parameter are coupled together, one of the two dynamical systems 
should be appropriately rescaled. Therefore, we propose a rescaled coupling scheme 

ciujw) (X {u — w), (2) 

where iD is a local average of w. It is important to understand that the coupling between two systems given by Eq. 
(0) is generally designed as a dissipative coupling. However, an interesting observation can be made at the limit of 
complete identical synchronization (i.e. the strong convergence of the two systems \\u{t) — w{t)\\ ^ for t ^ oo). 
From the point of view of image processing, the local average w is equivalent to the treatment of w by a low-pass filter. 
Moreover, at the limit of complete identical synchronization, {u — w) is equivalent to the treatment of w by a high-pass 



filter jlOJ. There is a similar effect on the second system under the same condition. As such, the rescaled sensor 
coupling expression given by Eq. (g) can be used for image processing, pattern recognition and shock capturing. In 
these applications, spatially selected treatment is of practical importance. To achieve spatial selectivity, we introduce 
following rescaled and adaptively distributed local sensors 

c(u, w) oc e{\ Ux \){u — w), (3) 

where the coupling strength e is a function of the gradient measurement | u^; |. For the purpose of edge-detected 
pattern recognition, we choose e(| u^ |) as a non-increasing function, e.g. 



£(| Ux I) = eexp 



2a2 



(4) 



where e and a are constants. For the purpose of noise reduction and oscillation suppression, we choose e(| Uj^ |) as a 
nondecreasing function, e.g. 

e(| Ux\) = e\ux |3, (5) 

where e is a constant. Obviously many other nondecreasing functions can also be used. In the rest of the Letter we 
restrict ourselves to the application of the present synchronization scheme to shock capturing. 

The solution of the inviscid Burgers' equation and the incompressible Navier-Stokes equations for very high Reynolds 
numbers is often difficult to attain due to the possible existence of shock front. Shock wave is a common phenomenon 
in nature, such as in aerodynamics and hydrodynamics, and is usually described by hyperbolic conservation laws 
and by inviscid hydrodynamic equations. The construction of numerical schemes that are capable of efficient shock 
capturing is a challenging task. 

To illustrate our synchronization approach for oscillation reduction, we first consider Burgers' model of turbulence 

Ut +UUx = -^Uxx, (6) 

He 

where u{x,t) is the dependent variable resembling the flow velocity and Re is the Reynolds number characterizing the 
size of the viscosity. Burgers' equation is an important model for the understanding of physical flows. The competition 
between the nonlinear advection and the viscous diffusion is controlled by the value of Re in Burgers' equation, and 
thus determines the behavior of the solution. We consider Eq. (o) using the following initial and boundary conditions 

u{x,0) = sm{TTx), M(0,i) = u(l,t) = 0. (7) 

The forth order Runge-Kutta scheme is used for the temporal discretization with a time increment At — 0.002. A 
discrete singular convolution (DSC) algorithm lOJl^ is utilized for spatial discretization with a total of 101 grid 
points in the computational domain. The use of DSC algorithm for solving differential equations has been extensively 
tested and further validation is given in TABLE I. 

Solving Burgers' equation at a high Reynolds number is a challenging task. At a Reynolds number of lO^^, the 
numerical solution quickly develops into a sharp shock front near x = 1. Severe oscillations occur near the shock 
front as shown in FIG. la. It should be pointed out that, almost all high order numerical schemes exhibit similar 
oscillations. To eliminate oscillations, we employ the rescaled adaptive coupling, Eq.(|5|), in solving Burgers' equation 
(O). Here w is computed by a local three-term average. Two systems, which are characterized by two Reynolds 
numbers (Rei = 1000 and Re2 = 100), are coupled with a coupling constant of e = —80. It can be seen from FIG. 
lb that all spurious oscillations could be eliminated. However, the synchronized solution is neither the true solution 
of Rc= 100 nor that of Re= 1000. Hence, it is desirable to have an oscillation-free solution at a given high Reynolds 
number. To this end, we design an auto synchronization approach by choosing two exactly identical systems, i.e., 
setting Rei =Re2 = co. ks two exactly identical systems are still coupled, oscillations are suppressed to a certain 
degree, depending on the coupling constant. For a relatively small coupling constant of e = —40, the solution is 
oscillatory at early times and become essentially non-oscillatory at a later time (see FIG. Ic). By increasing the 
coupling constant to e = —90, we have successfully eliminated all spurious oscillations as shown in FIG. Id. 

To validate the present approach further, we consider the two dimensional Navier-Stokes equation 

Uf + U • VU = -Vp + -^V^U + e{\ VU |)(U - W) (8) 

Re 

with the incompressible condition, V • U = 0. Here U — (u, w)"^ is the velocity vector, W is the velocity vector of the 
second system, p is the pressure, Reynolds number of Re= cxd defines the Euler equation. The domain of interest is 



a square [0, 27r] x [0, 27r] with periodic boundary conditions. Depending on the initial values, this system can be very 
challenging to solve. For smooth initial values, the solution scheme and the validity of the DSC algorithm were tested 
in Ref. ||I^. With appropriate initial values, the Euler equation (Re= cxo) can be used to describe the flow field of 
vertically perturbed horizontal shear layers around a jet. 

We now test our synchronization approach for the Euler equation with sharply varying initial values. This example 
is chosen to illustrate the ability of the present approach for providing very fine resolution even on a relatively coarse 
grid. The initial values are that of a jet in a doubly periodic geometry 



tanhl^^^ 



^ tanh(3^), ify>7r 
v{x,y,Q) = 5sm{x), (9) 



where 5 — 0.05 is used for the convenience of comparison with the previous study [ |l3yi4[| . This initial value describes 
the fiow field consisting of horizontal shear layers of finite thickness, perturbed by a small amplitude vertical velocity, 
making up the boundaries of the jet. However, this problem is not analytically solvable. Pioneering work was done 
by Bell et al ||I^ in this field, in which they utilized a second-order Godunov scheme in association with a projection 
approach for divergence-free velocity field with a general boundary condition. A state of the art high-order essentially 
non-oscillatory (ENO) scheme was constructed by E and Shu Um to resolve the fine vorticity structure of the double 
shear layers with periodic boundary conditions. 

We consider the parameter p — 7r/15, a case studied by Bell et al |13[ using a projection method with three sets 
of grids (128^, 256^ and 512^). E and Shu [|^ computed this case by using both spectral collocation code with 
512^ points and their high order ENO scheme with 64^ and 128^ points. The spectral collocation code produced an 
oscillatory solution at t = 10 (see FIG. 1 of Ref. ||ljl), while the high order ENO scherne produced a defect at i = 6 as 
the channels connecting the vorticity centers are slightly distorted (see FIG. 2 of Ref. |lj]). In the present simulation, 
we choose a 64^ grid for the computational domain with a time increment of 0.002. The synchronization prescription 
given in Eq. (||) is used with a coupling constant of e = —80. The results at different times (i=4,6,8, and 10) are 
plotted in FIG 2. It is seen that our solution is smooth (some non-smooth features in the contour plot are attributed 
to the coarse grid employed in the study) and stable for this case. In particular, no distortion is found in vorticity 
contours at t=6. For early times, present results compare extremely well with those of the spectral collocation code 
computed with 512^ points. There are no spurious numerical oscillations during the entire process. 

In conclusion, we propose the approach of synchronization as a robust, reliable and practical algorithm for shock 
wave computations. A length rescaling technique is introduced based on the dimensional argument of control process. 
To achieve computational efficiency and reliability, coupling sensors in spatially extended systems are adaptively 
selected according to local information. The coupling strength at each sensor is adaptively varied according to the 
magnitude of the local gradient of the system. The resulting coupled systems are analyzed from point of view of 
image filters. The proposed algorithm is validated by using Burgers' equation and the incompressible Navier-Stokes 
equation. A high accuracy discrete singular convolution algorithm [|ll| , [12| is utilized for the numerical simulation. 

For Burgers' equation, we have tested the computational accuracy and reliability at a moderately high Reynolds 
number. At high Reynolds numbers. Burgers' equation is difficult to solve. We first test the synchronization approach 
in two systems for different Reynolds numbers. One of the two systems develops spurious oscillations which could 
be eliminated by coupling to a low Reynolds number system. The oscillations are completely suppressed at synchro- 
nization. To make the algorithm practical for shock capturing at any given Reynolds number, we couple two truly 
identical systems (i.e., the same dynamical system, the same parameters). It was found that the oscillations could be 
completely eliminated above a suitable minimum coupling strength. 

To further validate the present synchronization approach for shock capturing in practical large scale computations, 
we consider the simulation of doubly periodic shear layer flows. For this purpose, the Navier-Stokes equation with 
appropriate initial and boundary conditions is employed. Numerical simulation of this system is an acid test for 
ordinary methods, including spectral methods. Our synchronization scheme performs extremely well. With a relatively 
small grid mesh, the present results are better than those of a state of the art shock capturing scheme, the high order 
essential non-oscillatory (ENO) scheme [^, obtained over a much larger grid mesh. This indicates that the proposed 
approach has a great potential for being used as an efficient and reliable algorithm for the practical simulation of fiuid 
fiows and computational physics in general. Progress has been made on the application of the present synchronization 
approach to image processing. 

This work was supported in part by the National University of Singapore. 
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TABLE I. I/i and Loo errors of the DSC solutions for Burgers' equation 



t 


0.4 


1.2 


2.0 


3.0 


10 


20 


40 


Li 


2.1(-4) 
2.8(-3) 


4.2(-5) 
6.8(-4) 


9.4(-7) 
1.2(-5) 


3.8(-8) 
4.0(-7) 


1.5(-11) 
4.0(-ll)) 


3.2(-12) 
5.3(-12) 


4.1(-13) 
6.4(-13) 



Figure Captions 

FIG. 1. Synchronization profiles of Burgers' equation at t =0.2 (i), 0.4 (ii), 0.6 (iii), 1.2 (iv) and 2.0 (v). (a) e = 0, 
solid line: Rei = 1000, dots: Rea^lOO; (b) e = -80, solid line: Rei = 1000, dots: Re2 = 100; (c) e = -40, Rei=Re2=oo; 
(d) e = -90, Rei=Re2=oo. 



FIG. 2. The vorticity contours of the synchronization solution of the 2D Euler equation with 64^ points. Up left: 
t = 4; up right: t ~ 6: low left: t — 8; low right: t = 10. 
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